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Abstract. We discuss the existence of inflationary solutions in a class of 
renormalization group improved polynomial f(R) theories, which have been studied 
recently in the context of the asymptotic safety scenario for quantum gravity. These 
theories seem to possess a nontrivial ultraviolet fixed point, where the dimensionful 
couplings scale according to their canonical dimensionality. Assuming that the cutoff is 
proportional to the Hubble parameter, we obtain modified Fricdmann equations which 
admit both power law and exponential solutions. We establish that for sufficiently high 
order polynomial the solutions are reliable, in the sense that considering still higher 
order polynomials is very unlikely to change the solution. 
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1. Introduction 

The idea that a perturbatively nonrenormalizable theory could be consistently defined 
in the UV limit at a nontrivial fixed point (FP), often called "asymptotic safety", is 
theoretically very attractive, especially when applied to quantum gravity [1]. Much 
progress in this direction has come from the direct application of Renormalization Group 
(RG) techniques to gravity. A particularly useful tool has been the non-perturbative 
Functional Renormalization Group Equation (FRGE) [2], defining an RG flow on a 
theory space which consists of all diffeomorphism invariant functionals of the metric 
g^v It defines a one parameter family of effective field theories with actions T^g^v) 
depending on a coarse graining scale (or "cutoff") k, and interpolating between a 
"bare" action (for k — > oo) and the ordinary effective action (for k — > 0). When 
applied to the Einstein-Hilbert action the FRGE yields beta functions [3, 4] which 
have made possible detailed investigations of the scaling behavior of Newtons's constant 
[10, 5, 6, 7, 8, 9, 11, 12, 13, 14, 15, 16]. It has been shown quite convincingly 
that the dimensionful Newton constant G is antiscreened at high energies, a behavior 
that eventually leads to the UV FP that is necessary for asymptotic safety. These 
analyses have then been enlarged to include matter [17] and a growing number of 
purely gravitational operators in the action. Truncations involving terms quadratic 
in curvature have been considered in [18, 19, 20, 21]; even higher operators could be 
included restricting oneself to powers of the Ricci scalar [22, 23, 24]. It has also been 
seen in the latter studies that the dimension of the basin of attraction of the FP (the 
so called "UV critical surface") is likely to be equal to three. Reviews of this work have 
appeared in [25]. 

The RG flow of the effective average action, obtained by different truncations of 
theory space, has been the basis of various investigations of "RG improved" black hole 
[26, 27] and cosmological [28, 29, 30, 31, 32, 33, 34] spacetimes. In particular, very 
recently it has been shown that the "RG improved" Einstein equations admit (power-law 
or exponential) inflationary solutions and that the running of the cosmological constant 
can account for the entire entropy of the present universe in the massless sector [35, 36] 
(see [37] for an extended review.) 

These works were based on the following logic [33]. If we want to study the quantum 
evolution of the cosmic scale factor, we should in principle use the full effective action 
r(g /MU ), which, as we mentioned above, coincides with T^g^v) in the limit k — > 0. 
However, our knowledge of this functional is rather poor. One way of gaining some 
traction on this issue is to observe that the Hubble parameter appears as a mass in 
propagators. Thus, the contributions of quantum fluctuations with wavelenghts greater 
than H^ 1 are suppressed. As a result, the functional at k comparable to H should 
be a reasonable approximation for the same functional at k — 0. We do not know Tk 
much better than the full effective action, but we can easily calculate the dependence 
of some terms in Tk on k. By doing so we effectively take into account nonlocal terms 
that would be very hard to calculate otherwise. This is analogous to what happens in 
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the discussion of the Coleman- Weinberg potential: there we restrict ourselves to quartic 
potentials but we identify the RG scale with the field itself. Since the quartic coupling 
runs logarithmically, this is equivalent to having a term 4 log(0 2 ) in the effective action. 

Previous investigations along these lines, in particular [35], have been based on the 
Einstein-Hilbert (EH) truncation. It is important to establish that the results obtained 
there persist when further operators are included in the truncation. We know that at the 
FP the coefficients of these terms are not very small, but their presence does not seem to 
affect the values of the cosmological constant and Newton's constant too much. In other 
words, the FP that is is seen in the EH truncation seems to be robust. The question then 
is to see if this stability of the FP against the inclusion of new terms is reflected in the 
stability of the corresponding solutions. This question is important because the values 
of the couplings at the FP are fixed and as a consequence there are no free parameters 
to be varied. We will see that the (power-law or exponential) inflationary solutions 
are indeed stable against the inclusion of new terms, but establishing this fact requires 
including a rather large number of terms. 

A somewhat different perspective on this subject has appeared in [38]. We will 
discuss in the conclusions the relation of that approach to the one adopted here. 

Our work is organized as follows. First we discuss the equations of motion in 
the presence of running couplings (section II), then we specialize the equations to 
Friedmann Robertson Walker (FRW) backgrounds (section III). Section IV is the core 
of the paper, where we explain in detail our procedure. The solutions are described in 
section V. Section VI defines a different "RG improvement" procedure and discusses 
the corresponding results. Section VII is devoted to a discussion of the apparent energy 
nonconservation in the "RG improved" theory. In section VIII we compare our approach 
to other work in the literature and we discuss the main aspects of this work that need 
further development. 

2. RG improved gravitational dynamics 

The effective action for the metric containing the effect of matter and graviton 
fluctuations is a complicated functional consisting of all possible integrals of scalar 
functions constructed with curvatures and covariant derivatives of curvatures. If the 
effective action is evaluated at some reference energy scale k, then the couplings 
appearing in this effective action will in general depend upon k. The dependence of 
the couplings on k is described by /3-functions, which collectively describe the RG flow 
of the theory. It is practically impossible to say something on the beta functions of all 
couplings. In this paper we will consider a subclass of terms that are somewhat simpler 
to study than the rest, namely functions of the curvature scalar R. This class of theories 
has been widely studied in the literature, at least at a classical level (see for example 
[39] for a detailed review). For the time being, in order to keep it as generic as possible, 
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we write the action in the form 



S = J d A x^\g~\F(R) (1) 

where F is some function of the scalar curvature. Note that even though we will refer 
to S as "the action", we always mean "the effective average action", meaning that 
fluctuations of the fields with momenta greater than k have already been integrated 
out. From (1) we obtain the equations of motion in the form 



where 



5S 



g\ Sg^ 



= F\R)R llv - ^F(R)g^ - V»V V F'(R) + g^V p V p F\R) (3) 

and we have written the energy-momentum tensor of matter on the right hand side. A 
prime stands for the derivative with respect to R. We want to use these equations to 
describe the cosmological evolution of the early universe. 

The "RG improved" cosmological equations will be obtained by replacing the 
gravitational couplings (cosmological constant, Newton's constant etc.) by running 
couplings. Insofar as the cutoff scale is identified with a function of the metric, and the 
metric is itself a function of the coordinates, the running couplings will also be functions 
of the coordinates. We will discuss some other general consequences of this assumption 
in section VII; for the time being we just proceed assuming that the function F appearing 
in (2) is built with such coordinate-dependent couplings, and that this dependence has 
to be taken into account when taking its derivatives. To be more explicit, if we assume 
that F(R) can be represented by a series of the form 

oo 

F(J2) = (4) 

i=0 

then 

V M F' = l)g t R^ 2 V^R + V^R 1 - 1 ] 

i 

While this seems to us the most natural and correct procedure, it is not a priori obvious 
that leaving out the last term (i.e. performing the "RG improvement" after having 
taken the derivatives) would necessarily be wrong. In order not to burden the reader 
with a doubling of all results, in most of this paper we shall follow the former procedure, 
and then devote section VI to a discussion of the results using the latter. Fortunately, 
we will see that several results are largely independent of this choice. 
In the following we will parametrize the action as follows: 

F = T^(f(R) ~ 2A) (5) 

107rG 
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Then (2) takes the following form 

f(R)R, u - \f{R)g^ - V,VJ'(R) + g^V 2 f'(R) + ^-V v f'{R) 

+ ^V»f'(R) - 2^V p f'(R)9„ - f'(R)G(v»V^ - g, u V 2 ^) (6) 
= - Ag^ 

We will set f{R) to be a polynomial of degree n, and write it as 

n 

i=l 

with fx = 1 by definition. Again, it is understood that when derivatives act on f(R), 
also the couplings fa have to be derived. 



3. RG improved Friedmann equations 

As we are interested in the cosmological evolution, we specialize to a spatially flat 
Friedmann- Robert son- Walker metric and take T^ v = diag(—p, p, p, p) to be the energy 
momentum tensor of an ideal fluid with equation of state p(p) = wp, where is 
a constant. Normally at very high energy it is natural to assume w = 1/3, and we will 
mostly do so. However, this is just an approximate description of the matter content of 
the Universe. As we shall discuss in section VII, due to the coarse graining the energy 
momentum tensor could have unusual properties and the effective w could be different 
from its classical value. 

In a FRW cosmology with scale factor a(t) we can write both G^ v and R^ v in terms 
of the Hubble rate H(t) = a(t)/a(t). In particular, we have 

Gtt = 3H 2 , R u — — 3 (H + H 2 ), R. u = R 

R = R» /i = Q(H + 2H 2 ), DR= -R-3HR (8) 

so that the (tt)-component and (minus) the trace of (6) become 

A{H) = 8vrGp + A (9) 
B(H) = MGp (1 - 3w) + 4A (10) 



where 



A{H) = - 3{H + H 2 )f + 3Hf> + If - 3H^f (11) 




B(H) = -6(H + 2H 2 )f + 2f + 3f'+\ 

f , GG - 2G 2 + 3HGG 
61 G* {U) 
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One can eliminate p from (10), thus obtaining an equation that determines a(t), 
while (9) is used to determine p: 

B{H) = (1 - 3w)A(H) + 3(1 + w)A (13) 

The equations for the case without matter can be obtained by setting p = in 
(9,10). Then a(t) can be obtained by solving 

4A(H) = B{H) (15) 

while either (9) or (10) provide an additional equation that involves A. As we shall see, 
this system is quite constraining. 



4. Cosmology in the fixed point regime 

4-1- Fixed point action 

The RG flow for F(R) theories of gravity has been studied in [22, 24, 23]. One can 
actually derive a beta functional for the entire function F, but the corresponding FP 
equation is very complicated, so the analysis has been done by expanding F in Taylor 
series. It has proved possible to study truncations involving up to eight powers of R. 
We now recall the results of this analysis, in the parametrization provided by equations 
(5) and (7). 

The asymptotic safety scenario posits that in the ultraviolet the couplings reach 
a fixed point. This statement has the obvious meaning when applied to dimensionless 
couplings such as the electromagnetic coupling, or the quartic coupling in scalar field 
theory. In the case of dimensionful couplings, it means that they must tend to constant 
values when measured in units of the cutoff. We call k the cutoff, and we denote by a 
tilde the ratio of any quantity by the cutoff raised to the canonical dimension of that 
quantity. Thus quantities with a tilde are by definition dimensionless. For example 
R = R/k 2 is the curvature measured in cutoff units, and 

A = A/k 2 ; G = Gk 2 - f % = f t k 2l ~ 2 (16) 

are the couplings measured in cutoff units. In the following we will refer to them as 
"the dimensionless couplings". It is paramount to understand that the existence of 
a fixed point refers not to the dimensionful couplings such as A, G etc. but to their 
dimensionless counterparts A, G etc. As a result, at a fixed point the dimensionful 
couplings do run quite violently, namely they depend on the cutoff precisely as dictated 
by naive dimensional analysis. This is to be contrasted to the low energy regime we are 
familiar with, where the opposite happens: the dimensionful couplings do not run, and 
consequently the dimensionless ones have a very strong classical running. 

We now come to the fixed point of F(R) gravity. The values of the dimensionless 
couplings at the fixed point, for various polynomial truncations of order up to ten, are 
listed in Table 1. The first eight lines are taken from [24]. For reasons that will be 
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0.9583 
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7 


0.1202 


0.9488 


1.61 


466.6 


501.8 


288.4 


-163.1 


-218.6 








8 


0.1221 


0.9589 


-4.225 


413.6 


430.2 


328.1 


-56.17 


-298.6 


-226.3 






9 


0.1242 


0.9715 


-13.17 


398.5 


362.9 


297.7 


28.77 


-223.2 


-219.3 


-127.7 




10 


0.1242 


0.9718 


-13.16 


391.1 


360.7 


307.8 


3.874 


-230.1 


-229.3 


-123.9 


40.32 



Table 1. Position of the FP as a function of n, the order of the truncation. To avoid 
writing too many decimals, the values of /* have been multiplied by 1000. 



discussed in section V.B, these truncations are actually insufficient for our purposes and 
we have found it necessary to examine also the truncations n = 9, 10. All the other 
details of the calculation (cutoff, gauge etc.) are the same as in [24]. We have found 
only the two fixed points listed in the last two lines in the table. We have not calculated 
the critical exponents that pertain to these fixed points, so our understanding of their 
nature is less complete than for the truncations n < 8, but the close resemblance of the 
values of most couplings is a strong hint that these are indeed the correct prolongations 
of the n = 8 fixed point to higher truncations $. 

In this paper we want to use this information in cosmology. Since the universe 
becomes hotter and more strongly curved as one proceeds backwards towards the big 
bang, it is clear that in order to explain the dynamics of the early phases of its evolution, 
we need to have a theoretical model for the behavior of interactions at high energies. 
Asymptotically safe gravity is a model for what happens near the Planck scale: it 
posits that all interactions reach the fixed point regime. So, assuming that gravity is 
asymptotically safe, in this paper we want to study the evolution of the cosmic scale 
factor in the very early universe, when the (dimensionless) couplings are so close to their 
fixed point value that we can actually use the values give in table I. 

4-2. Cutoff choice 

The RG equations give us the dependence of the couplings on some cutoff scale k, but 
they do not tell us what k is. In order to apply them to cosmology, or any other problem, 

I One may be worried by the sign flip of /g. In this connection we note that the coefficients arise 
from the sum of a large number of terms and that there is nothing to guarantee their sign. The actual 
difference between /g in the truncations n = 8 and n = 9 is 0.08; similar shifts occur also elsewhere in 
the table. 
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one has to identify a physical quantity that acts like a cutoff. From a Wilsonian point 
of view, the cutoff is a typical energy or momentum scale of the phenomenon under 
study, and one integrates out all fluctuations of the fields with momenta higher than 
the cutoff. But this still falls short of identifying it uniquely. The best choice of k has 
to be made on a case by case basis, analyzing the physical problem at hand. 

In our case, as also pointed out in [38], the Hubble constant H provides a natural 
infrared cutoff in the loop diagrams of the gravitons (see also [37] for an extended 
discussion of this point). We will therefore assume that the cutoff is 

k = £H (17) 

where £ is a positive number of order unity S. 

The logic that we will use is then as follows: in the fixed point regime we can replace 
k by £H in (16). This turns the dimensionful couplings A, G etc into functions of time. 
We then use these expressions in the dynamical equations (9) and (10). The dimensionful 
couplings are replaced by powers of H, dimensionless couplings which we take from table 
I and the free parameter £. The occurrence of the time dependent function H, where 
previously there appeared a constant, obviously complicates the equations significantly. 
The resulting equations should give a reasonably good approximate description of the 
dynamics of the universe in the fixed point regime. We then look for de Sitter type 
solutions 

a(t) = aoe Ht ; H = constant , (18) 
or power law solutions 

a(t) = a t p ; H = V - . (19) 



4-3. Reliability analysis 

One important advantage of having the results of increasingly more complete truncations 
is that this gives us some quantitative handle on the reliability of the calculations. To 
understand this point let us write 

n 

f(R) = k 2 f(R) ; f{R) = Y,lR l 

i=l 

The fixed point polynomial f*(R), whose coefficients are given in table I, is plotted in 
figure 1 for n = 2,4,6,8,10. Clearly the approximation becomes more accurate as n 
increases. It is worth emphasizing that table I does not give the Taylor expansion of a 
fixed function: if that was the case, then a single row of coefficients would have been 
enough. Instead, for each n, the truncated fixed point equations give a whole new set 

S Different choices of cutoff have been considered in the past, for example k ~ l/a(i) [40, 41] and 
k ~ 1/t [28, 33, 34]. The choice made here has been considered previously in [35]. See also [42] and 
[43] and more recently [44, 45]. Note that for a power-law dependence of the scale factor the choices 
k ^ 1/t and k ~ H arc equivalent. 




of coefficients, which provide a polynomial approximation for the "true" fixed point 
function /*. The fact that the numbers in the columns of table I do not change too 
wildly is an encouraging sign that the data in the table do indeed resemble a Taylor 
expansion of some function. As expected, one sees that the shape of the function near 
the origin is rather unchanging, and that the uncertainty moves progressively to larger 
R. In order to make this quantitatively more precise we shall use the following method. 
Let /* be the fixed point curve in the truncation n. We say that /* is reliable as long as 
it differs from f* +1 by less than Af*. Here we will take Af* = 0.025: since /* ~ 1 for 
R S3 1, this seems a reasonable criterion for two successive approximations being "near" 
to each other in some domain. Then we find that the truncations are reliable for 

R < c (20) 

where c is given by the following table: 



n 


1 


2 


3 


4 


5 


6 


7 


8 


9 


10 


c 


0.45 


0.32 


0.45 


0.87 


1.03 


0.86 


0.76 


0.87 


0.99 


0.99 



This method cannot give a reliability value for the highest truncation, so we 
conservatively assume that the n — 10 solution has the same range as the n = 9 
one. 

Using (8) and (17), equation (20) implies that 

ce>12 + 6^- 2 (21) 

Later on, when we discuss solutions, we will make sure that they occur within these 
bounds. 

One may wonder how things change when one uses different definitions for the 
"nearness" of two functions. Of course if one chooses Af* to be too small, then also the 
best truncations become unreliable and conversely if Af* is too large then all truncations 
are reliable. Values of Af* between 0.01 and 0.1 give similar results. The values of c in 
the table have a clear tendency to increase, but they have also some randomness. This 
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can be decreased by choosing a different criterion for what one means by nearness, for 
example defining c by the requirement that the integral of the modulus of the difference 
of two successive approximations be less than a preset value. We think that the simple 
criterion used above is sufficient for our purposes. 



5. Cosmological solutions 

5.1. Einstein- Hilbert truncation 

In order to make contact with [35] we discuss first the case n = 1, corresponding to the 
Hilbert action f(R) = R. In this case the expressions (11) and (12) reduce to 

A(H) = 3H 2 - 3H^ (22) 

Gr 

B(H) = 6(H + 2H 2 ) - ^G-2G 2 + 3HGG ^ 

Gr 

It is not difficult to show that power law solutions of the type (19) are obtained for 



V J 1-H, (24) 

3-A^ 2 

In particular we see that p is real and positive for £ greater than a critical value. 

In reference [35] a different RG improvement was used: the running couplings were 
inserted in Einstein's equations, in such a way that the G terms were absent. We defer 
a detailed discussion of the difference between these procedures to section VI. For the 
time being we merely notice that using the approach described in [35], for w — 1/3 one 
obtains 

p = i (25) 

6-2A4 2 V ! 

The exponents are plotted in Fig.2 as functions of £. Note that the positive branch 
of (24) is qualitatively similar to (25) (the short-dashed curve). Unfortunately both 
curves are in the region where the truncation is not reliable according to the criterion 
(20) (they are in the grey area). It is interesting to notice that for £ > 4 the negative 
branch of (24) is within the domain of validity of our approximation, but then the 
exponent is too small for inflation. 

5.2. The case n = 2 

As another explicit example let us give the form of the equations when / = R + /^-R 2 - 
Then we have 

A(H) = 3H 2 + 18f 2 (2HH - H 2 + QH 2 H) 

+ 3Qf 2 H(H + 2H 2 ) - 3H^ \l + l2f 2 (H + 2H 2 )] (26) 

Gr L 

B(H) = 6(H + 2H 2 ) + 36/ 2 (# (3) + 7HH + 12HH + AH 2 ) 
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Figure 2. Black solid curve: Exponents p of power law solutions as a function of £ for 
n = 1, w = 1/3; see equation (24). The short-dashed curve is the solution (25), found 
in [35] and further discussed in section VI. The long-dashed curve indicates p = 1. 
The grey region represents the region excluded by the reliability criterion, see (21) 
with c = 0.42. 



+ 36f 2 (2H + UHH + 6H 3 ) + 36 f 2 (H + 2E 1 



- 72 



G 



f 2 (H + 2H 2 ) + f 2 (H + 4HH) 



-3 
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GG - 2G 2 + 3HGG 
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l + 12f 2 (H + 2H 2 



(27) 



The power-law exponents p are now obtained from the solution of a cubic equation. 
Their explicit form reads 



where 



and 



36 



Pi 



P2,3 



1 + w) (l8f~; + e) + 2 -j^ + 2 2/3 M 2 /3 



18(1 + w)^ 2 (3- A*£ 2 
72(1 + w) (l8/| + e) ' 24/Hl ^f 3)Ml - 2 2 / 3 (l T < V3)-Mf 
18(1 + w)£ 2 



Mi = 162(1 + iy) ^2(1 + w) ^18/ 2 * + £ 2 

-e 2 ("3(11 + 3^)/! + e 2 ) h-Le 



M 2 = 11664(1 + w) 3 (18/* + «e 2 J - 8748(1 + w) 2 £ 2 ( 18/* + £ 2 
x (3(11 + 3^)/* + £ 2 ) (3-A4 2 

52488(1 + w) 2 fZ^ (3 - ie 2 " ' 



(29) 
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Figure 3. Black solid curve: Exponents p of power law solutions as a function of £ for 
n = 2, w = 1/3; see equation (28). The dashed curve indicates p = 1. The grey region 
represents the region excluded by the reliability criterion, see (21) with c = 0.24. 



+ 



-AM\ + 8503856(1 + wf 4(1 + w) ( 18/* + £ 



(30) 



- 3£ 2 (18/* + e) (3(11 + 3^)/* + e) (3 - A*£ 

1/2 



18/*£ 4 (3-A^ 



The exponents as functions of £ are then plotted in figure 3. It turns out that p\ is 
negative and therefore uninteresting for cosmology; The other two solutions are complex 
in general, but for £ greater than some critical value they are real. We see that these 
solutions are very similar to the ones found in EH truncation. Of the two branches, 
only the lower one is reliable for large £, but again it is too small to be of interest for 
inflation. The upper branch is not in the region where the truncation is reliable. 



5.3. The general case: de Sitter solutions 

If we make the ansatz (18), the terms containing /' and G in equations (11,12) vanish, 
and one finds 

B = 4A . 

Inserting in equations (9,10) there follows that w = —1, i.e. matter must have the same 
equation of state as the cosmo logical constant. We may then as well absorb such matter 
in the definition of A and set p = 0. Thus without loss of generality we will study the 
exponential solution only in the absence of matter. 

Since R = 12H is constant, equation (13) reduces to 

i?/'-2/ + 4A = . (31) 

This equation had been studied earlier in [46]. We render the equation dimensionless 
by going to tilde variables and use the values of the dimensionless couplings given in 
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table I. Then we solve for R. Numerically we find solutions for n — 2, 3, 6, 7, 8, 9, 10. 
These are given in the second column of table II. Note that for any solution R*, the 
value of the parameter £ is fixed. This is because in de Sitter space R = 12H 2 , so 
R* = R*k 2 = R*H 2 !; 2 implies that £ 2 = 12/ R*. The solutions mentioned above would 
correspond to the values of £ shown in the third column in table II. 

No real solutions are found in the truncations n = 4,5. In order to understand the 
reason for this consider the coefficients in table I and observe figure 4 where we have 
plotted the left hand side of the equation. The behavior of the function Rf — 2f for 
large R is determined by the sign of the coefficient (n — 2)f*. For n = 2, this leading 
term cancels and the equation is determined by the linear term which is negative, so 
that there is a zero. For n > 2 the sign of this coefficient is the same as the sign of /*. 
For n — 3, it is positive and the curve is bent upwards but not enough to eliminate the 
solution. For n = 4, 5, fl and fg are again positive and the curve is bent upwards enough 
that the solution disappears. If this was as far as one could get with the truncation, then 
one might conclude that the solution that is present for n = 2, 3 is a truncation artifact. 
In fact the solution occurs at values of R that are outside the reliable range defined by 
(20). Now, when we add the terms of order 6, 7 and 8, the coefficients of the highest 
terms are negative and therefore a solution reappears, first at a rather large (and hence 
unreliable) value of R but then at reasonably small R. Is this sufficient evidence for the 
existence of the solution? If we were to rely only on the results of [24] , where n = 8 
was the highest truncation considered, one would feel that the evidence is somewhat 
inconclusive. In fact, looking at figure 4, where n = 8 is represented by the dashed curve, 
one can easily imagine that if fg was sufficiently positive the solution could disappear 
again. It is for this reason that we have looked at the fixed point condition in the cases 
n — 9, 10. Luckily fg is again negative, and /* is positive but not very large, so that 
a solution exists in all these cases and actually occurs within the domain of reliability 
of the truncations. Furthermore, the position of the solution seems to be quite stable 
for n = 8, 9, 10, which suggests that the higher order terms will not affect it too much. 
The hard lesson that one learns from this is that it may be necessary to go to very high 
truncations before one obtains reliable results. 

5.4. The general case: power law solutions 

Now we look for power law solutions (19). In this case the condition (21) implies 



When one uses (11), (12), (7), (8) and (17), for dimensional reasons these equations 
reduce to the product of some power of t times a function of p and £. Solving them for 
all t means that p and £ must satisfy two algebraic equations, and one expects to find 
at most isolated solutions. At least for n = 2 no solutions were found. 




(32) 



Consider first the case without matter. Then we have to solve the equations 
A(H) = A 3B(H) = AA 



(33) 
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Table 2. Table 2. De Sitter solutions for various truncations. When solutions exist, 
only the smaller one is displayed; no solutions are found for n = 4, 5. 




Figure 5. Solid curve: numerical solutions for the exponents of power law solutions 
as a function of£ for n = 10, w = 1/3. The dashed curve indicates p = 1. The grey 
region represents the region excluded by the reliability criterion, see (21) with c = 1.06. 



Inflationary solutions in asymptotically safe f(R) theories 



15 



Let us now search for power law solutions in the presence of matter. Things now 
look better because matter gives us a new degree of freedom and allows us to find 
solutions for continuous ranges of values of £. Namely, for fixed £ one can use equation 
(13) to determine the exponent p and then use (14) to fix p(t). In this way one avoids 
having to fix £. The case n = 2 has already been discussed in section 5.2. The cases with 
n > 2 cannot be solved analytically, but solutions exist and can be found numerically. 
As before we use w = 1/3. 

Doing this highlights once again the need to go beyond the truncation n = 8. For 
suppose that we only knew the results up to n = 8. Then according to our criteria we 
would have to take for this truncation the same reliability range as the n = 7 truncation, 
which has a smaller value c = 0.79. With this criterion, the whole n = 8 solution that 
asymptotes to de Sitter would be unreliable. It is only by going to higher n that we 
can validate the solution in an acceptable range. We note that part of these curves lie 
still to the left of the reliability limit (32), but that is the part that is less interesting 
physically. As a confirmation of the results, we show in figure 5 the exponents p in the 
truncation n = 10. 

6. A more restrictive RG improvement 

Let us return to equation (2). As we mentioned, there is a possible ambiguity concerning 
the stage at which one should replace the usual "constant" couplings by time dependent 
"running" couplings. To illustrate this point, let us consider for a moment the special 
case of the Einstein-Hilbert action, F = j^q(R — 2A). If we were to simply "RG 
improve" Einstein's equations, as in [35], we would have 



where G is now allowed to depend on x. This is different from the procedure that we 
followed in this paper, which leads to the general equations 



For want of a better terminology we will refer to these two procedures as the "restricted" 
or the "extensive" improvement respectively. At least at the level of the Einstein-Hilbert 
truncation it is not a priori clear that one of these procedures is right and the other 
wrong. For example, one may derive Einstein's equations from arguments unrelated to 
an action, then the first choice seems a legitimate one. 

Equation (34) has been written before in [33], where the restricted improvement 
is referred to as "improving equations" and the extensive improvement as "improving 
actions". This terminology is unambiguous in the case of the Einstein-Hilbert action, 
but it is ambiguous when one considers more general actions such as (1): in fact is quite 
clear that also the extensive improvement could very well be described as "improving 
equations", when the equation is written in the form (2). Also, while our extensive 





(34) 
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improvement may well be seen as replacing the constant couplings by position-dependent 
couplings in the action, the procedure that we have followed has been quite different 
from the one advocated in [33], because we have not treated the position-dependent 
couplings as prescribed external functions. Having hopefully clarified the differences in 
these various approaches, in the absence of a very strong a priori argument against the 
restricted improvement, in this section we will describe the results that one obtains from 
it. The equations of motion are of the form (2) where, now, 

Ep, = F'{R)R^ - ^F(R)g^ - F"(R) (V M V,i? - g^V 2 R) 

- F"\R) (V^RV V R - g^VRf) (35) 

Let us now discuss solutions. First we observe that in the search of de Sitter 
solutions, the difference between the extensive and the restricted improvement is 
immaterial. To see this compare the difference between the definition of B in the two 
cases. This difference is given by 

- it q fGG ~ 2G 2 + 3HGG 
, G / ~ 6f CP 
where do means that one only takes the derivative with respect to t of R and not of the 
couplings. Each term in this expression contains at least one time derivative of H and 
therefore vanishes for de Sitter space. Thus, also with the restricted improvement the 
de Sitter solutions are given by table 2. 

When we look for power law solutions, additional solutions appear for n > 2. We 
consider first the case n = 2. 

A(H) = 3H 2 + 18f 2 (2HH - H 2 + QH 2 H) 

13(H) = 6(H + 2H 2 ) + 36/ 2 (^ (3) + 7HH + 12HH + AH 2 ) (37) 
One finds that there exist power law solutions with exponents 

4 (54(1 + w)f* + £ 2 ) + + 2^(M> 2 ) V3 

P'l = 7 — 2 — \ (3* 

6(1 +w)£ 2 (3-A^ 2 ; 



3d 2 f + 9Hd f - 6^f - 3f ^ ™^r»»™ m 



P2,3 



8 (54(1 + W )f* 2 + e) ~ 2i/Hl ^f 3)Ml ~ 2 2 '\l T i V3)Ml /3 



where 



and 



12(1 + w )e (3- 

M[ = 4 (54(1 + w)ft + - 54(1 + w)(ll + 3w)f;e (3 - A*?) (39) 

M' 2 = 16 (54(1 + w)f* + - 324(1 + w)(U + 3w)j;e 
x (54(l + ^)/* + e 2 ) (3-A4 2 



+ 1944(1 + w) 2 f* 2 ^ (3 - A4 2 



2 
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Figure 6. Solid curves: numerical solutions for the exponent p of power law solutions 
as a function of £ for n = 10, w = 1/3, using the restricted improvement. The dashed 
curve indicates p = 1 . The grey region represents the region excluded by the reliability 
criterion, see (21) with c = 1.06. 



+ 



-A(M[) 6 + 16 4 54(1 + + £ 



(40) 



81(1 + w)(U + 3w)f* 2 e ( 54(1 + w)f* 2 + e ) ( 3 - A,£ 

1/2 



-486(1 + w)7 2 T 3-A4 



The remarks in the end of section V.B apply here too. Solutions (38) are all real 
for sufficiently large £; p\ is always negative, while p 2 diverges for finite £ and P3 stays 
almost constant and is too small for inflation. 

For higher truncations the solutions have again to be found numerically. Figures 6 
gives the exponents as functions of £ in the cases n = 8 and n — 10 respectively. We 
note that there is a solution that starts at p ~ 4 for £ = and has the same de Sitter 
asymptote as the one found with the extensive improvement. These solutions are very 
close in the domain of reliability of the truncation. There is then another solution which 
starts at p ~ 1/2 for £ = and has p > 1 for 3.15 < £ < 4.75, and whose physical 
meaning is doubtful. We conclude that the physically most relevant part of the solution 
is rather insensitive to the choice between extensive and restricted improvement. 



7. Energy and entropy 

One somewhat unsettling aspect of this approach to cosmology is non-conservation of 
the energy momentum. From equation (2) one finds that 

2V fi E» v = V fl T\ (41) 

and both sides of the equation would vanish if they were obtained by varying a 
diffeomorphism invariant action. But the RG improved equations were not simply 
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obtained by varying a diffeomorphism invariant action: after the variation, the couplings, 
which are usually treated as constants, were replaced by functions of the metric. Of 
course if we replaced the couplings by scalar functions of the metric in the action, 
before varying, then we would obtain another diffeomorphism invariant action. But this 
is not the procedure that we use here, so one should not expect the l.h.s. of (41) to 
be zero. As a consequence, also the r.h.s. cannot be zero, so the energy momentum 
tensor on the r.h.s. of the RG improved equations cannot be obtained from varying 
some diffeomorphism invariant matter action. 

Let us calculate the l.h.s. of (41). Since is linear in F, it is easiest to do this 
when the function F has a Taylor expansion as in (4). Then one finds 

V»E» V = ~%F (42) 
where V means that the derivative acts only on the couplings and not on R: 



oo _ , oo 

v uF = v„ gi R i = ^ E ( 43 ) 

i=0 i=0 

This equation has a very simple interpretation: the failure of energy-momentum 
conservation (a gravitational anomaly) is proportional to the beta functions of the 
couplings. 

At this point it is important to stress that we are not claiming to have found 
a violation of energy momentum conservation at a fundamental level. If, as we said 
in the introduction, we were using the full effective action (namely the functional Yj. 
at k — 0), then there would be no RG improvement and, barring the occurrence of 
genuine gravitational anomalies (as discussed for example in [47]), there would be no 
anomaly because the full effective action is diffeomorphism invariant. So the type of 
gravitational anomaly that we are discussing here is entirely due to the RG improvement 
and the associated coarse graining. 

To further discuss this point we can recast (6) in the form 

R w ,-\gp,R = T™ + T lu , (44) 



where 



rpKG _ X 

f'(R) 



\gM{R) - Rf'(R)) + V M V„/'(i2) - g»vV 2 f{R) 



V,GV u f'{R) V u GV,f{R) WGV p f{R) 
G f'(R) G f'(R) G f'(R) 9 ' w 1 ] 

A 



f'(R) 



and T M j, = (8ttG/ f'(R))T llv . It is thus clear that a non-vanishing Tj^p exactly 
compensates the nonconservation of T^, so that the total effective stress-energy tensor 
T*°* = T^, + T^y given from the r.h.s of (44) is conserved, in general. 

One can think of T^ u as the stress-energy tensor for a dissipative, i.e. non-ideal 
fluid which interacts with the coarse grained gravitational degrees of freedom. In fact in 
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the functional modes with wavelengths smaller that 2%/k are integrated over. Since 
k is identified with H, it decreases with time and therefore, as time proceeds, more and 
more gravitational modes are being removed from the description of the system. Such 
modes carry energy and it should not be surprising that when they are removed energy 
seems not to be conserved. 

Although in this paper we have considered the coarse-graining of the modes of 
the gravitational field, similar considerations would apply also to matter fields. It 
is interesting to review how this happens already in the more familiar framework of 
Wilson's action for the Ayj 4 /4! scalar field theory. Quite generally, the Wilsonian action 
in this case can be defined as 

e"^* 1 = J D[ip] ]j6(Mx) ~ $(x))e~ 5M = J D[ip]e- S " lv ' 9] (46) 

where <fik{x) is the average of the field <p in a domain of characteristic length 1/k. For 
actual calculation it is possible to introduce a smearing function Vk(x, x') which is nearly 
constant within distances shorter than 1/k but rapidly decays to zero outside this region, 
so that we write 

<f> k (x) = / d d x u k (x,x')ip(x') (47) 



One can evaluate the on-shell condition for the blocked action (46) by means of the 
standard saddle-point approximation = ^Jx 1 ^ (see [48] for details) which gives 



Dtp = ^(p 3 + 2M 2 ( J d d x' u k (x, x')if(x') - $ 



(4£ 



The non-local contribution coming from the coarse-graining kernel acts as a source terms 
in the equation of motion for the effective theory at the scale k. This is responsible for 
a modification of the standard conservation law which now reads 

V„T/ = 2M 2 9^(0 fc - $) (49) 

where T^ u is the stress-energy tensor of the original (bare) field ip. It should be noticed 
that in the limit k — > 0, (j) k — $ and one recovers the standard conservation law with 
no additional source term. 

Let us now see in detail how the modified conservation law arise in the case of the 
Einstein-Hilbert truncation. The RG improved Friedmann equations read 

3H 2 - = 8vrGp + A (50) 

, • On GG — ~h 3UGG / v / h. 

6(H + 2H 2 )-3 — = 87rGp(l-3w;) + 4A (51) 

Gr 

Taking the derivative of the first equation and using again both equations, one is led to 
the following modified continuity equation 

p + 3H(p + p) = V (52) 
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where 



V 



8nG 



A + 8npG 



3G 
8ttG 2 



H 



HG 
~G~ 



H 2 



(53) 



The l.h.s. is just V^T^q, so this equation agrees exactly with the time component of 
(41), when we use the definitions 

2A 1 



go 



9i 



16%G ai 16ttG 

Is is important to stress that this expression can consistently be obtained from the 
4-divergence of Eq.(44), as it must be for consistence. Note that we do not use the 
additional consistency condition described in [33] as we allow for an unobstructed 
exchange of energy between matter and the gravitational effective dynamics. The 
first term in square brackets in the r.h.s of (53) had already appeared in [35], but 
the second one is due to the different improvement scheme used here. We will discuss 
these differences in some detail the next section. Here we recall from [35] that the time 
variation of the couplings, and in particular of the cosmological constant, can be seen 
as a transfer of energy to the matter degrees of freedom and hence gives rise to an 
increase in entropy of the cosmological fluid. In fact it was shown that essentially all 
of the entropy that is observed can be accounted for in this way. We note that this 
calculation, which is given in equations (3.9)-(3.15) in [35], does not depend on the 
form of V as a function of A and G, so the entropy generation will be the same in the 
approach followed here. In fact, as discussed in [35], the peak of the entropy production 
is in the crossover region, near the Gaussian fixed point, where only the running of the 
cosmological constant is relevant and G ~ 0. 

We now give the formula for energy (non) conservation in the general case. In 
general we can write 



oo 



4HE, 



00 



HE* 



Then using 



oo 



A(H) - A 
IQttG 



B{H) - 4A 
IQixG 



and the explicit expressions (11,12), one finds 



V = 2V„E% 



8nG 



A-A§-i(/-/'*) + i/§ 



This agrees with the result of inserting (5) into (42). 

We can see the implications of these facts for the search of cosmological solutions 
at the fixed point. At a fixed point we have 



so inserting in (42) we find 



A = (4 - 2i) 9i 



£(4 - 2%^ = -2^(2F - RF>) 



i=0 
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There are therefore two ways in which the anomaly could vanish: the first is that 
2F = RF', which is satisfied if and only if F — g2R 2 - This is because #2 is dimensionless, 
and therefore it is constant at a fixed point. The other way, which could work for any 
form of the function F, is that V u k = 0. This depends on the choice of the cutoff, 
i.e. how it depends on the metric, and on the solution. In particular, if we look for 
vacuum solutions T^ v = 0, we must also have V^T^j, = 0. Vacuum solutions must have 
a vanishing anomaly and this is why they are harder to come by than solutions with 
matter ||. If we make the cutoff identification k = £H, and if F is not just g2R 2 , then 
a necessary condition to have a vacuum solution is H = 0, so the only vacuum solution 
is de Sitter. 

The (non) conservation of the energy momentum tensor with a restricted 
improvement works in a different way. Instead of (42) we find now 

= R^F' - Ul v F - (V^V U R - g, u V 2 R)WF" 

- (VpRVvR - g^(VR) 2 )V^F"'(R) (54) 

Note that only the second term is present if we perform the extensive improvement. The 
new terms that are seen here can be viewed as additional contributions to the anomaly. 

Finally let us discuss the energy density that is required for the existence of the 
power law solutions. Once the solution for a(t) is found, one can plug it in (14) to obtain 
p. It turns out that p always depends on time as t~ 4 , as demanded by dimensional 
analysis, so the only issue is the value of the constant prefactor. Explicit formulae are 
easily derived in the case n — 1. One finds 
_ A{H) - A 
P ~ 8ttG 

3w + A^ 2 ± J 3(1 + w) (s(w - 1) + 2A*£ 5 



9^ 



(55) 



in the extended RG improvement scheme and 

= 3# 2 -a = sie 2 i f56) 

in the restricted improvement scheme. Note that the additional term proportional to 
G, which is present in the former case, is negative, making the density p negative in the 
former case and positive in the latter. 

These facts signal that a detailed analysis of the matter sector is necessary. We 
shall not undertake this analysis here, leaving it for future work. We just add a few 
remarks. It is not difficult to see that one can have positive energy also with the extended 
improvement, provided w is smaller than a critical value w cr < —1. Furthermore, in 
a more refined treatment we should perform the coarse-graining also on the matter, 

|| The same argument holds more generally if we demand that the energy momentum tensor be derivable 
from a diffcomorphism invariant matter action. 
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treated as quantum fields. On the other hand, it can be argued that the realistic RG 
trajectory is the one which emanates from the non-Gaussian FP and spends a long time 
near the Gaussian fixed point, corresponding to the classical era [35]. During this time 
G ~ and one recovers the more familiar "restricted" RG evolution where the density 
is always positive. 

8. Discussion 

The results obtained here can be summarized by saying that if we identify the cutoff 
with a multiple of the Hubble parameter, as in (17), inflationary power law solutions 
(i.e. with exponent p > 1) exist for 2.7 < £ < 3.9. The dependence on £ is strong, with 
the exponent diverging at £ ~ 3.9. Exactly at that point, the theory admits a de Sitter 
solution. Provided that £ lies in the above range, that the starting point is close enough 
to the FP and that p > 1, it should always be possible to have a sufficient number of 
e-foldings. 

Weinberg [38] has studied the FRW equations that follow from a general 
gravitational action containing arbitrary powers of curvature, and the possibility that 
they admit inflationary (de Sitter) solutions. In his approach the cutoff is a fixed mass 
scale that has to be optimized for the treatment of inflation. Unlike the approach 
followed here, it does not depend on time and therefore there are no "RG improvement" 
terms in the equations. This may sound like a very different procedure, but in practice it 
is not, for two reasons. The first is that the fixed optimal cutoff is tuned to a description 
of inflation, and if one wanted a cutoff that is tuned to some later stage of the cosmic 
history, it would be different; thus effectively one would have again a time-dependent 
cutoff. Conversely, when we focus only on de Sitter solutions, also our time-dependent 
cutoff k = £H becomes time-independent. Therefore, if one started with the equations 
of [38] and truncated the action to the form F(R) that we consider here, then one would 
find the same solutions. 

The equations of motion determine if as a function of the cutoff, so if we assume 
that the cutoff is a multiple of H, the coefficient £ is determined by the solution, as we 
have observed. It is worth noting that the value of £ that produces a de Sitter solution 
is in the right physical range, namely the cutoff retains fluctuations with wavelengths 
that are few times smaller that the horizon scale. 

In this work we have discussed the effect of performing the RG improvement in the 
equations of motion, after having varied the action. To see how things could work in a 
different approach, let us momentarily view equation (31) as a differential equation for 
/ (instead of an algebraic equation for R, as we have looked at it so far). It has the 
solution A = 0, / = g2R 2 , which corresponds to a scale invariant theory. (We have also 
seen in section VII that in this case the anomaly automatically vanishes.) For such an 
/, de Sitter space is a solution for any value of R, and £ 2 = 12/ R remains undetermined. 
We observe that this theory can be obtained if we make the RG improvement in the 
action, before varying. In fact, if we make the cutoff identification k 2 = uR, in the fixed 
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point regime ^ = 2% R 2 \ so 

OO 00 

i=0 i=0 

Note that since the fixed point coefficients c/i seem to be all less than one, and to have 
alternating signs, with u> > 1 there is reasonable hope that the series converges. The 
same behavior obtains with the cutoff identification (17), if one restricts onself to de 
Sitter spaces. This result also agrees with the expectation that the fixed point should be 
described by a scale invariant action. Cosmological solutions of this theory have been 
studied in [49]. 

There are properties of the theory that are not universal, but depend on the choice 
of the cutoff. Clearly, it will be important to check that the observable low-energy 
physics is not strongly dependent on the regulator choice, but this question can only be 
addressed with a proper numerical integration of the complete system of improved RG 
equations and the /3-functions for the couplings. 

One aspect of the problem that we have not touched upon here is the end of 
inflation. In [38] this is estimated by considering the development of instabilities in the 
exponential solution. In the approach that we followed here one would assume that the 
initial point of the RG evolution is not at the fixed point but somewhere close to it. The 
solutions described here would then hold only approximately. The RG evolution would 
take the theory away from the fixed point, with a speed that can be calculated from 
knowledge of the scaling exponents of the theory. Eventually when one gets sufficiently 
far from the fixed point the inflationary solution would give way to an ordinary radiation 
dominated universe. One would think that by choosing the initial point sufficiently close 
to the fixed point one can have an inflationary period of arbitrarily long duration. Of 
course the two approaches should give at least similar results. We hope to return to this 
point in the future. 

Unfortunately the results of this paper are not yet a realistic basis for a model of 
inflation. The main reason is that we are neglecting a great number of terms in the 
action. For the study of Robert son- Walker cosmologies, Weyl terms are unimportant, 
and for n = 2 (meaning with four derivatives) there is nothing beyond R 2 that is not 
a total derivative. However from n = 3 upwards there are many terms in the action 
that contain traces of powers of the Ricci tensor, whose effect we are currently unable 
to estimate. One rather sobering result of our analysis has been that in order to find 
reliable solutions one has to go to really high orders in the derivative expansion. It is 
encouraging, however, that the (technically unreliable) results of the Einstein-Hilbert 
truncation [35] proved in the end to give the correct qualitative picture (at least within 
the class of F(R) truncations). One may hope that also when Ricci terms are included, 
the results of the low order truncations are not too misleading. 
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